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Most often, an odd prime number has its "half size" 
number, (P-1)/2, composite. For some primes (e.g., 47), 
the half size number is also prime, and the quarter size 
number, (P-3)/4, may also be prime. The pattern on the 
cover shows how this fact can generate an interesting path. 


The rule of formation is: if the half size number 
is composite, advance in a straight line. If the half 
size number is prime, but the quarter size number is 
composite, turn to the right; if both the half size and 
quarter size number are prime, turn to the left. 


This path must eventually cross itself. 


The Problem is: What cells will contain more than 
one prime, and what are those primes? 
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Think Before You Cross San Pasqual 


---by John Todd 


California Institute of Technology 
Special Freshmen Lecture, January 27, 1965 


The activities of the Numerical Analyst are in some 
ways similar to those of the Highway Patrol but the analogy 
is not perfect and should not be driven too far. The 
Numerical Analyst tries to prevent computational catastro- 
phes by ensuring that reasonable and prudent procedures 
are used and by defining a safe computing code. These 
are in the nature of service functions: as a researcher 
he is on the lookout for the exploitation of computers in 
new areas. 


It is the business of the Computing Center to see 
that good equipment and advice is available and that 
prices are within our reach. 


It is certainly not the purpose of the Highway Patrol 
to prevent people using the roads; similarly, if there was 
no computing going on, I would be probably doing what to 
me would be less interesting mathematics. 


I shall today try to show you by simple examples, 
most of which can be done by paper and pencil, some of the 
pitfalls of computation. A little imagination will 
enable one to guess what can happen in current practice, 
where millions of operations take place in a typical 
calculation. 


I hope these examples will frighten you, but not 
too much. The scientist who disregards the computer today 
is foolish and puts himself at a disadvantage with his 
competitors. Equally, the scientist who blindly uses 
the computer (and this means the accompanying software more 
than the hardware) is asking for trouble. 


i. We now learn about associativity in high school and 


1 


that abc” ~ can be evaluated either as 


first ab, then (ab)cn2 


1 


or as first be, then a(be7). 


¥* othe undergraduate houses are on the south side of San Pasqual and 
our computer center is on the north side.” 
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If we imagine the operations of multiplication and 
division as those used by the computer, with rounding to 
a fixed precision (for example, to two decimals) this is 
no longer the case. If we take 

asec be .1l c= .13 


then the two calculations are: 


.0132 rounded to .01 divided by .13 gives .076 
which rounds to .08; and 


.11 divided by .13 gives .846 which rounds to .85 
which multiplied by .12 gives .104 which rounds to .10. 


The correct answer is .10. 


2. Consider the solution of the system of equations: 
| (et SxS + 6x3 + 5x, = 23 
8x, + 6x, + 10x, + 9X, = 33 


| 7X4 + 5X5 + 9X3 + 10x), = Sple 
By inspection the exact solution is 

x, = 1, x5 21, x3 21, xy 21. 
It is also easy to verify that 

X1 = 9.2, Xo = -12.6, X35 4.5, xy = -1.1 

nearly solves the system, indeed up to errors of .1, -.1, 
-1l, -.1 in the right hand sides. Indeed, if we perturb 
the right hand sides to 

324 €, 23 -€, 33+ €, J1 = € 


we can verify that the exact solution to this new problem 
is 


X, 5 1+ 82e 
Xp ® 1 - 136€ 


1+ 35e 
Xy =il1- 2le 


Relative error in the data is magnified by several 


’ thousand in the solution. 


(Due to T. S. Wilson) 


2.2 (Due to W. Kahan) 


A similar phenomenon occurs in the 
case of the system: 


-2161x + .1441y - .1440 
1.2969x + .8648y - .8642 


ut 
oO 


6) 


The exact solution of this system is 


Be Wo oS 
However, 
x = .9911, y = -.4870 


satisfies these equations up to errors of 1078 Sores 


é} 


Sic In paragraph 2 we saw the effect of small changes in 
linear equation problems. We now consider the same 


effect in polynomial equations. 


(3.1) Bez Ge li eo =O 

has four roots: 1, 1, 1, l. 

If we change the middle coefficient from 6 to 
6 = 49 x 10° 


(that is, by less than 1 in Le the roots change by 
about 3 in 100, being: 


1.02681, 0.97389, 0.99965 + 10.026455 


(3.2) Similarly, if we change from the equation 


PLP non Cour a Ones 


-1 
the roots change from O to numbers of modulus 10 . 
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This can be interpreted in matrix language as follows: Let 

A be the 10 x 10 matrix with 1's in the (i, 1+1) positions 

fe OTN en A1VeN egies ogi) This matrix has the characteristic 6 
polynomial ,10 and all its characteristic values are 0. 

By changing only the (10,1) element from 0 to 10719, the 
characteristic polynomial of the perturbed matrix becomes 


tO - 10710 and the characteristic values become 


-1 ra rr 
10 cos — + i sin— 
( 5 5 ); 
® S Oy s5009e (Due to G. E, Forsythe) 
(3.3) The above examples suggest, and rightly, that 
multiple roots are especially sensitive. However, if 
we take an equation with roots 1, 2, ..., 20, say 
20 
Il (z-r) = me Ne he eae (20!) = 0 
rel ® 
19 -23 
and change the coefficient of z from 210 to 210 + 2 
(that is, we make a relative error of about lone then 


the roots of the new equation are 
1.00000 0000 
2.00000 0000 
10.09526 6145 + 0.64340 Og044 


3.00000 0000 
11.79363 3881 + 1.65232 97281 


4.00000 0000 
4.99999 9928 
6.00000 6944 
6.99969 7234 
8.00726 7603 
8.91725 O2k9 
20.84690 8101 ” 


13.99235 8137 + 2.51883 00701 
16.73073 7466 + 2.81262 48941 
19.50243 9700 + 1.94033 03474 


(Due to J. H. Wilkinson. This example was also discussed in issue 23 
of POPULAR COMPUTING .) 


4, Again, consider calculating sin(7/6) on our 2D 
computer. 


We have to replace 7/6 by .53, replace sin x by 
x - (1/6)x3 and obtain after several roundings: 


Sintt/6) = .5e 1.00 = .17( 52 x .52) 


S! (2 te = oes 
= 49, 
There are three types of error here. We could 


distinguish two other types before these: a wrong 
formulation of the problem: 


y' + y = O instead of y" - y = 0 


and an analytic error in choosing sine instead of a 
cosine. 


Se Newton's method for obtaining square roots is 
motivated as follows. if x, is an approximation to N 


say, too low (high), then N/xXp will be an approximation 


which is too high (low). Hence their average 


al 
*nt1 aun $ (W/7n)} 
should be a better approximation. 


For instance with N = .25, XQ = 1 we obtain 


X, = 0.625, Xp = 5.125, x3 = 0.500152, x, = 0.500000... 


It can be proved that if €, = x, -VN then 
2 -1 
x 
n 


aoe 
S Be emis i 


n+1 


--the error at one stage is of the order of the square at 
the preceding stage, which means that, roughly speaking, 
we double the number of correct decimals in our answer 

at each application. 


It can also be proved that xp4j -x, = 3(N - se) <0 
provided that XQ >0, n>O so that Xn41 < %, and that x 


does converge to WN. See the diagram. 
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5.2 (W. G. Hwang, J. Todd) 


The recurrence relation 


Mery os xi oN = x2)/2N 


which also converges to the square root of N for suitable 
Xg Was popular in the days of computers without division 
since it does not involve division by a variable quantity, 
as does the relation just discussed. This relation 
converges almost as fast as the earlier one. Take N = 2, 
Then with a good guess we have 


Xo = 1.5000 
2a ad 1.4062 
Xp = 1.4141 
x3 = 1.4142. 
However, this relation has curious properties. It is not 


unexpected that a bad guess leads to trouble. With a 
starting value of 3.000, we find successively 


X, = -2.2500 


X5 = -0.5273 
X3 = -0.7543 
Xe = -1.4142 


so that we have convergence to -/2 instead of 2 ! 
However, a worse guess gets to the right answer: 


the starting value 3.1000 leads to the values -2.7978, 
+1.2782, +1.3952, +1.4138,...,+1.414213. 


A very bad guess leads to great trouble: x, = 4 leads to 
values of -10, 235, -3244116, 8.5355 x 10 t6 the 18th power. 


The investigation of the behavior of this sequence 


when x, is taken in the range V6 <x < Vi0 is very instructive: 
infinitesimal changes in XQ can change the limit from 


+/2 to -v2 . 
[see also the formula given in PC22-2: 
Xnt1 = -5%y(3 - Nxf) 


which involves no division at ail. | 


Quadratic convergence toVN 


All this 1s theoretical arithmetic; when we come to 
practical computation the infinite descent guaranteed by the 
strictness of the last inequality just cannot happen. The 
question arises: at what stage do we stop and take the current 
Xy, as the required square root? It can be shown that the 
appropriate x, is that one which first satisfies Xn 2 Xyn-13 


that is, when the sequence becomes stationary, or reverses 


its direction. Let us look at two examples, the first on 
our 2D machine. 


(1) N= .01, X%9 = .11, N/xg = .09 
x, = -50X, + -50(N/X) 
= .06 + .05 
re os Xo 


This example shows that strict inequality need not happen. 
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(2) N= -1/2, = Al 


a6) 
ie 245 2-1/9, -17/112. 


Hence, ¥-1/2 = -17/112. 


This example shows the necessity to check that N20, 


Naturally, if one is asked directly to fing wee » one takes 
appropriate action; in practice, however, intermediate steps of 
calculation are rarely monitored by humans, and the subroutines 
must have safeguards incorporated. 


There are now many books discussing, at various levels, 
numerical mathematics, both theoretical arithmetic and 
practical computation. There are also various collections 
of algorithms (in books and in periodicals) and various 
implementations of these are generally available. 


It is always advisable, before beginning any serious 
computation, to carry out controlled computational experiments 
with the programs contemplated; by this we mean to do test 
problems whose exact result is known and observe the accuracy 
obtained. For this purpose, collections of test problems 
in various areas are available. 


Although it is not generally possible to give realistic 
error estimates for non-trivial problems, it is often possible 
to combine the results of computational experiments and the 
technical knowledge of the customer to make error statements 
which have some authority, but not that of traditional 
mathematics. 


Although it is fun to compute, it is often quicker 
and cheaper to get the result from tables and all scientists 
should be familiar with such books as 


A. Fletcher, J. C. P. Miller, L. Rosenhead, and 


Addison-Wesley. 


National Bureau of Standards, Handbook of 
Mathematical Functions, edited by M. Abramowitz 
and I. A. Stegun, U. S. Government Printing Office. 


(At $11.50, the latter book is undoubtedly the best bargain 
in books available to the scientist today, with 1045 large 
pages of mostly useful material.) 


@ 
_] 
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The square root recursion: 


x = x(3N = x©)/2N 


ntl 
in Professor Todd's paper (section 5.2) has interesting 


behavior, as he points out. The behavior can be 
summarized this way: 


0 ve Q Vio 


Diverging 
0 O 0 oscillating 

oO 

[ee] 

i cea 

= 

The regions above the horizontal line are values for xo, | 

the starting values to use in the recursion. Below re] 

the horizontal line are the values to which the process S 
goes, either converging (for starting values less than the o. 


square root of 10), or oscillating (at the square root of 
10), or diverging (greater than the square root of 10). 


The number Q bears looking into. It is approximately 
3.037. Just what is this value? 
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6 
be 


FREE FRPWWWW WWWWWWWw 


.7084.2976926618947 3 
«73251115681724824 3 
-75628575422107 2007 
.77976 3149684619494 
.802952460761391619 
825862 365544778203 
.848501131276805069 
.870876640627796747 
. 89299641587 3260547 


-914867641168863596 
.9 3649718310217 3196 
-957891609680405479 
.979057207896 391860 
.000000000000000000 
.020725758589057976 
.041240020622190271 
.061548 100445679789 
.081655101917 348071 
.101565929702347522 


4 .121285299808556820 
4 .14081774942285 3250 
4 ,160167646103808229 
Horepeeete oso ueen ee 
te 
i 
hk, 
4, 


198336453808407722 


.217163326508746214 
-23582358425489 3164 


4 ,254320865115005776 
272658681697916825 
29084.0427026207112 


The N-series, a feature of POPULAR COMPUTING 
since issue number 2, will be discontinued. 
Although it has been favorably received, its 
usefulness, particularly for the less common 


functions, has attenuated for most readers. 


FREE EFF 


FREE E FE EEL 


. 308869 380063767444 
. 3267487 10922225147 
- 344481485768611902 
. 36207067 1454837565 
- 379519139887889266 
. 396829672158179276 
-414004962442103773 
.43104762169 3634159 
-447960181138631042 
-4647450955845 376 34 


-481404746557164709 
~497941445275414797 
~514357435474001380 
.530654896083492777 
- 54683594 377634 3894 
. 5629026 35386966728 
.578856970213327471 
.594700892207039806 
.6104 36292058446570 
.626065009182741793 
.641588833612778893 


CUBE ROOTS 


The 


values for square and cube root, and the common 


and natural logarithms, are given here for 


’N = 51 to 100. 


.141428428542849997999 399811 3672652787 6617 11599027 
.2111025509279785862 384425 3494099189250259 31476905 
.28010988928051827 1097 30249 15270327937776696825765 
. 348469228 3495 34294591852224117674175897 84 24419700 
.4161984870956629487 113974408007 1 3060979904 3190975 
~48331477 35478827 711674974646 3309860 35120396155575 
. 5498 344. 35270749697 2 3668480694611705822219470462 34 
.61577 310586 39082856614 11027 158323005 3607055925466 
.68114574786860817 5769687 0217 3137 247 30624510776149 
- 7459666924 148 337 703585 3079956479922166584 34105832 
.810249675906654 3941297227 3575910141 3568 3051366486 
-8740078740118110196850 3444881200786 36810861220209 
-93725393319377177 15048472609 17 781277130777 5492474 
-00000D0000000000000000000000000000000000000000000 
062257 748298549652 36661 32 303037711 31134396 3056086 
. 1240384046 35960 36045988 35682660403485042040867253 
.18535277187 244996995 37037247 339294588804868154980 
.24621125123532109964.2819711948154050294 3984507472 
. 3066238629 18074852584.2627449074920102 322142489557 
. 366600265 3407554797817 2025785187489 392815 36929867 
.42614977 31763586 306 341399062027 360 316080024015608 
.485281374238570292810132 34525818847 14180 312522617 
-54400374531753116787 1648 3262 397064 34594455 3295333 
.60232526704.262677172947 35 350497 1 36 32027 5355572907 
660254037844 3864676 37231707529 3618 347 140262690519 


-717797887081 347 10447 396 3967 7 19231 3182740078504649 
-774964 387 392122060406 388 307416 30956087 587 6827 5545 
. 831760866 3278468547 640427 2695925 39641746 394809 314 
888194417 31558885009 14416754087 278170764506037295 
-9442719099991587856 3669467492510494.176247 34 384461 
.0000000000000000000000000000000000000000000000000 
.0553851 3813741662657 38081669840664130521244640969 
.110433579144298881945626 1046886691900991 391682650 
#165191 3899116800131 760913871 5601697 719509131 pana 
.229544457292887 31000227428176279 315724680504 

27 36184954957037525164 1607 399017462626 34689120763 
. 327 37905308881504555447 554232055698 3276240694.1917 
. 3808 315196468591091 312602270889 325611764567068235 
.43398113205660 381132066037 76226407 1698 36226334151 
.4868329805051379959966806 332981556011586654179757 
-5393920141694564915262158602322654025462342525055 
.59166 30466254 3908 3194876128 325 38783999 341408 38083 
.64.36507609929549957600310474 3266 318 390690 369 306 33 
PEE Gee SAE ETRE NEE Gore ie eee 


7224 


6794 34480896 39068 384131998996002992525839003375 


-1979589711327 12392789 13629882 356556786 37899226267 
-8488578017961047217462114149176244816961 362874428 
-8994949 366116653416118210694678865499877031276386 
-949874 37 1066199547 3447982100120600517812656367681 


SQUARE ROOTS 
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COMMON LOGARITHMS 


1.70757017 60979 363658 35197797 5834523392077 
1.7160033436 347991596 33982947 391314484 3661 
1.7242758696007890456 3299229 16272565926955 
1.73239 3759822968507098822604.4898 38954 3686 
1.740362689494 24 384 55 364610765185312149385 
1.748188027006200416 3534 3294276611527 37881 
1.75587485567 2491 3988 31 3613790120446271513 
1.76342799 35629 37 28254.658565769 37480180225 
1.770852011642144190260656 3845351442389267 


1.77815125038364 3632508766797979608 3359683 
1.7853298350107670338857485 13757 321 3492634 
1.7923916894982538748804429948429087490719 
1.799 34054945358170530227 20651028668118838 
1 .80617997 3983887 1712824 333683469581606091 
1.81291335664285557 399276626 32178 354040615 
1.81954 393554186867 32589667692226 325776750 
1.8260748027008264 341491 316292260685809496 
1.8325089127062 3631896764768 3777 32 30835439 
1.8388490907 37255316162805015506 3048588976 


1.8450980400142568307122162585926 3619 34836 
1.8512583487190752860928294 350354291 352704 
1.857 3324964 3126846023127 24906837096987048 
1.863322860120455901074 3869004703085 344529 
18692317197 30976192022189584.26 36224747512 
1.8750612633917000468675501138061292556637 
1.88081359228079 135196 38112652059153714875 
1 .8864907 2517248187 14624162298 356604 351903 
1.8920946026904804017 1527 19559219 367667980 
1.897627091290441427994821 3864782496864829 


1.90308998699194 3585641216684.17 34790803046 
1.9084850188786497491801116130204612368005 
1.91381385238371668972315074469267 38262987 
1.91907809237607 3903832760352027 2612470016 
1.9242792860618816584 347219512967 375562201 
1.9294189257142927 333264 309996038440032394 
1.9344984512435677216188270479537151855770 
1.939519252618618524627874666224 3703004544 
1.94448267215016862639 14 166554165033220113 
1.949 39000664491278472354 33697024411246652 


1.9542425094 39 3248745900558065102306184003 
1.95904139232109 35999187 214165 34964624 3133 
eae a ole ot oeopac ena > 4901 (001 | 20022358 
1.96 482911895 32392169017 3200337 39310315038 
1.973127853599698659627958294 17 369 36669280 
1.977723605288847 766 3225945810324 362911829 
1.9822712330395684 133637 2237687758044 30411 
1.9867717 342662448517 84 3618116655774494258 
Aol eacoT peek 620s 81714119097654137 353 
1.9956351945975499153402557777532548601070 


Bho 
.95124 37185814 27 3548879516844 8167174095626821 347100 
-97029191 35521218 34 1444691 390290577 70359977752901 34 
988984046564 274 38 36029678 322257 5 368201797180781850 
.007 3331852 324709186627029111913169 39 347 30820819578 
-0253516907 351492333570491078177094 3386 35851 32652 32 
0430512678 34550151404 27 2668810 37924 188486981911146 
06044 3010546419 336600504 15 382008817 357001 304827246 
-07753744 39057 1945061605037 37 1969762406 3346789 31976 


- O94 3445622221006848 304688 1 3065066480 32409218080094 
- 11087 3864.17 331124875138910 3425614746 31568174 307026 
-127134385045091555 346 39644.60005337785254 3906482974 
- 1431347263915 3268789584 3217 2882311389 3206584521612 
- 15888 308335967 18565033927 2874905940845 300080615016 
. 174 3872698956 37 110654246774791506244 330869299017 24 
.18965474.20264255448744.209 36 34 58 3157254469754610882 
- 204692619 39096605967007 1996 36 37 227 505669 3290 321014 
-21950770517610669908 39988607 894796717 39203281294 34 
- 234106504597 259 3822019980687 327 2182308987087 265114 


3 
3 
3 


east 


FO SAR aaa Or he Gt STC Ree OSE ad ase FEHR FEEL FERRE 


NATURAL LOGARITHMS 


9318256 327 24 3257 71644779854 7956522402 3569 357039892 


- 24849524 2049 35898912 334419812754 39 372 3818621819844 
- 26267987 704 1 3154.21 329454 5 3251 303409675957652669824 
-276666119016055 3110421868 3821958111352148151871378 
- 290459441148 3911290921088574 385425709047 5284485898 
- 30406509 3204169753785327 79248962 37 3197557772151514 
«317488113536 31044059676 390337490098 369869 3266 34686 
- 3307 33340286 33107884 3491674804 20667 3388 3795 3000660 


- 34380542185 368 3849 167296 321408 30902945879158 350612 
- 3567088266895917 36865964 799946020877 5282586 3692994 
« 369447852467021494 17 294554148141092217 354122440922 


. 3820266 3467 388161226968781905889 391182760189169602 


- 3944491546724 387655809809476901028185899622 31277 34 
406719247 26425 311328 3995494495584 15645191060374578 
- 418840607 79659792 347547 222 329137045 3029 31 3056649 32 
.43081679884 331 361533506222 32820585704 3557555611110 
-442651256490 31645485029 3951099 314175113804 36684012 
«45434729625 35077 328900746 348040236036 346 3631922754 
~ 465908118654 58 37 18578517 2692844 37 3101420034717 3108 
477 33681447820647 231363994233965900404820725700 368 
4886 36 3697 321 398 3831781554.0669849219404660 387118 32 


-4998096703302650668084819285294 1561689608260425950 
-51085950651685004 11588401850084983 34442352674 32718 
-521788577049040 30964 12170747 26549254 59 338058 354 596 


. 53259949 31532559 37 32440956 14648829 150974294 8828824 
.54.329478227000 38962 381827912 3035027697 155063636524 
. 55387689 16005408 3460978676511404117676298061555228 
. 564 3481914678 362 3848140584421 340854502499122960850 


- 57471097850 33828221 1672162170396171 380891490264 312 


.5849674786705719 196279 37608 3445 36027 34966959 350818 
.§9511985013458992685 Ban 32091810180 7031 .CE eels 
-60517018598809 13680359 2909 3687 284.15202202977 24154 
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Problem solution 


Contest 14, Generating Triangles (PC45-1), called 
for analyzing all combinations of 3 points taken on a 
10 x 10 grid. With the coordinates of the 3 points 
taken as 2-digit numbers, each combination produces a 
4th point by summing the squares, modulo 100. For example, 
the combination 2,2; 4,8; and 8,3 leads to: 


22% = 48h 
48° = 2304 
83° = 6889 


9677 = 77 mod 100. 


The 161,700 possible combinations were to be analyzed and 
tabulated as follows: 


a) A triangle is formed for which the fourth 
point lies inside the triangle. 


b) A triangle is formed for which the fourth 
point lies outside the triangle. 


ec) A triangle is formed for which the fourth e 
point lies on the triangle. 


ad) The three points do not form a triangle (or, 
the triangle has zero area). 


Winner Andrew Vettel, Jr., West Mifflin, Pennsylvania, 
wrote his program in BASIC for a HP-9830 calculator. After 
a 50 hour run, he obtained the results: 


a 9187 
b) 139179 
e 8886 
d 4448 


Mr. Vettel's approach to the problem was to generate 

the 4th point and calculate the areas of the three triangles 
formed between that point and the original three, taken two 

at a time. If one of these triangles has zero area, then 

the 4th point lies on the original triangle. If the sum 

of twice the areas of the three new triangles exceeds twice 

the area of the original triangle, the 4th point is exterior; 
otherwise, it is interior. J 


Problem 170 (PC49-20) presented a sequence that 
began as shown: 


Three questions were posed: 


1. What is the algorithm for 
producing the sequence? 


2. What is the 1000th term? 


3. What fraction of the terms 
are perfect squares? 


The first two questions have been answered. 
Dr. N. J. A. Sloane, of the Bell Laboratories, 


reports that his colleague, C. L. Mallows, 
found that: 


7 t. 
ten a E ioe . 1)| 


where the square brackets denote largest 


integer. This recursion can be expressed 
as: 
t4, =+(3t, - 2) if t is even 
nt1 2 n n 
a 1 
tit =o (3t, - 1) if t, 1s odd. 


Meanwhile, Prof, James L. Boettiler, 
Talladega College, Talladega, Alabama, 
utilized a high precision package he had 
written, and calculated the terms of the 
sequence to the 1000th term: 


13344 15470 41507 96868 59187 

72545 97117 77188 63580 57690 

80462 19427 45785 80740 25513 

67719 51737 06881 64437 02825 

33207 82974 37156 87939 78924 

31324 07622 85342 14982 05714 

28739 41052 43478 73072 31331 (ea 
93 
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136216568 
204 324851 
306487276 
459730913 
689596369 
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MATRIX INVERSION 


A symmetric matrix can be inverted by the following 
procedure. The given matrix is A, which has been 
partitioned into four sub-matrices. The inverse matrix 
is C, similarly partitioned. The procedure calls for 
inverting the upper left sub-matrix and one other matrix 
of the same size. 


e Alo) he Cia” 
LAs, Ago- Cai So2. 
ES 
M= Ajq Ajo 
=} 
Coo = [Aap - Ap, Mi] 


Cees eae 
aed 

Cay > “a0 
= le 

Crier slit 21 


If the given matrix is non-symmetric, multiply it 
on the right by its transpose (thus forming a symmetric 
matrix). Invert that matrix, and multiply the result on 
the left by the transpose of the original matrix. All 
of this may be expressed as: 


ms ie 
Airs AGAMA a 


The product of the original matrix and its inverse 
is then: 


ew ee io can Ayi°q2 + 42°27 


ei 
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Ao2%a1 Ao1°10 + Aaa@aa - 


which will be the identity matrix, as required. 


For testing purposes, a matrix can be constructed 
by taking a matrix of the form: 


1 6) 6) 0) 
A 1 6) 0 
B C 1 fe) 
D &: F 1 


and multiplying it on the right by its transpose. The 
result will be a symmetric matrix; moreover, if the 


variable elements are integers, then the resulting matrix 
will have an inverse with integer elements. 


For example, a test matrix is constructed according 
to the given plan: 


TG 10. 0.) [¥>-2 5 ; 
Ce OF Olle lo 2.4 6 2 

Pet weolmlo) Gx 1 74. _alg 
eee Wi u|O 0 02 5 | 


and the steps in inverting that matrix are these 


= : ey 


ero) <2 
M «= = 
ae 5 ee 1 
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_ fe = @ 
i iy 

Fis oa | 20% — »<28 

-7 aL }-158 22 

(| 

2 

Bo ari fool, —=15e) 814 -638 

Riese). | 28. 122 638 501 ® 


814 -638 201 -28 
-638 501 -158 22 
201 38-158 50 -7 
-28 22 -7 1 


